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1  Introduction 


Various  physical  phenomena  are  observed  to  occur  in  the  arctic  and  its  subregions.  The 
most  noted  one  if.  the  sea  ice  itself  that  reflects  a  phase  change  of  sea  water  in  response  to 
atmospheric  cooling.  Once  formed  it  acts  as  a  insulator  between  the  ocean  and  the  ice. 
The  resultant  oceanic  structure  beneath  the  ice  acquires  a  distinct  characteristic  relative 
to  the  open  water  regions.  Temperatures  beneath  the  ice  tend  to  be  near  freezing,  while 
the  salinity  and  with  it  the  mixed  layer  dynamics  are  controlled  by  the  freezing  and 
melting  of  the  ice  through  the  subsequent  brine  injection  or  dilution.  The  mixed  layer 
depths  themselves  tend  to  vary  from  30  to  70  meters  under  the  ice. 

The  arctic  ocean  can  be  compartmentalized  naturally  into  different  regions,  depend¬ 
ing  on  the  predominant  processes.  In  the  upper  layer,  turbulent  mixing  and  Coriolis 
effects  are  equally  important  and  as  a  result,  this  layer  has  Ekman-like  dynamics.  In  the 
interior  the  pressure  gradients  tend  to  balance  the  Coriolis  force,  resulting  in  a  quasi- 
geostrophic  flow.  At  the  bottom,  there  is  a  benthic  boundary  layer  that  has  Ekman  type 
dynamics  and  is  forced  by  the  interior  flow.  A  lateral  transport  of  vorticity  proceeds 
from  the  basin  boundaries  toward  the  interior. 

In  such  an  approach  the  resulting  oceanic  velocity  can  be  thought  of  as  a  superposi¬ 
tion  of  velocities  arising  from  the  various  regions:  top,  interior,  and  bottom.  We  exploit 
this  property  by  essentially  solving  separate  models  for  the  various  regions  of  the  ocean, 
with  coupling  between  them. 

The  upper  ocean  is  modelled  in  the  framework  of  a  three-dimensional  mixed  layer 
approximation,  and  is  coupled  to  the  Hibler  thermodynamic-  dynamic  ice  model.  Two 
different  modelling  approaches  are  used  for  the  interior  ocean.  In  one,  the  geostrophic 
velocity  is  obtained  from  an  inverse,  ^-spiral  type  model.  In  another,  the  barotropic 
velocity  is  calculated  prognostically.  Topography  is  included  in  both  approaches. 

The  seasonal  and  regional  behaviors  of  the  ice-ocean  structure  is  studied,  along  with 
the  temperature  and  salinity  distributions  (under  the  ice,  in  the  MIZ,  and  in  the  open 
water)  In  addition,  the  nature  of  the  mixing  processes  and  their  effect  on  the  upper 
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ocean  heat  and  salt-  budgets  are  anlyzed.  Finally,  we  compare  the  ice  edge  locations 
against  analysed  results  derived  from  remote  sensing  and  some  in-situ  reports. 

The  seasonal  behavior  of  the  barotropic  flow  has  also  been  modelled  in  the  Arctic, 
including  the  Greenland  and  Norwegian  Sea  basins.  The  time  behavior  of  the  associated 
transports  has  also  been  studied. 


2  Modeling  Approach 

For  ice  representation  the  dynamic-thermodynamic  model  of  Hibler  (1979,  1980)  was 
used.  In  this  model  the  momentum  balance  of  the  ice  is  expressed  as 

m(Du/dt)  =  mf(k  x  v)  +  ra  +  fh  -  mgV H  +  F  (1) 

where^ m  is  the  ice  mass  per  unit  area,  u  is  the  ice  velocity,  /  is  the  Coriolis  parameter, 
ra  and  fj  are  the  atmospheric  and  bottom  stresses,  H  is  the  sea  surface  dynamic  height, 
and  F  is  the  force  due  to  variation  in  the  internal  ice  stress. 

The  ice  thickness  and  compactness  equations  are  treated  as  in  Hibler  (1979),  together 
with  the  Semtner  (1976)  thermodynamics  model  for  ice  growth  and  decay. 

The  equations  describing  the  upper  mixed  region  are  as  follows 

tiT/tit  =  0I$z{~1FT'  +  K(t?T/dz))  +  (l/pU)c)(i?F/i?z)  -  V  •  (HT)  +  AV7T)  (2) 

OS/xH  =  d/i?z(-^  +  *(t?S/t?z))  -  V  •  (uS)  +  A  •  VJS)  (3) 

x)u/x)t  =  d/t Jz(  —  w'u'  +  k(t9u/\ )z))  +  fv  —  Du  (4) 

dvjdt  =  \ ?/tfz(-tt;'v'  +  K[dv/dz))  —  fu  —  Dv  (5) 

where  T  is  the  mean  temperature,  5  is  the  mean  salinity,  u,  v  and  w  the  x-,  y-,  and 
z-components  of  an  advective  oceanic  current  velocity  (to  be  defined  later),  F  the  down¬ 
ward  flux  of  solar  radiation,  pv  a  reference  water  density,  c  the  specific  heat  of  seawater, 
D  a  damping  coefficient,  k  an  ambient  diffusion  coefficient,  /  the  Coriolis  parameter,  A 
the  horizontal  eddy  diffusion  coefficient,  t  the  time,  and  z  the  vertical  coordinate  (posi¬ 
tive  upward  from  the  level  sea  surface).  Primes  indicate  departure  from  the  means.  The 
advective  terms  are  retained  in  the  temperature  and  salinity  equations,  but  neglected 
in  the  momentum  equations.,  based  on  a  scale  analysis. 

The  advection  current  is  given  by 

u  =  uj  +  u„u  =  v,-  -f  vt,  w  =  W{  (6) 

where  the  u,-  and  vt  are  components  of  the  instantaneous  wind-drift  current,  u>,-  is  the 
vertical  motion  resulting  from  the  divergence  of  the  U(  field,  and  ti,  and  v ,  are  compo¬ 
nents  of  a  geostrophic  current. 

The  Level-2  turbulence  closure  method  of  Mellor  and  Yamada  (1974)  is  used  to 
parametrize  the  vertical,  eddy  fluxes  of  temperature,  salinity  and  momentum.  The 
geostrophic  velocity  is  calculated  with  a  /?-spiral  type  inverse  method  of  Peggion  (1988). 
The  model  assumes  the  usual  hydrostatic  and  geostrophic  approximations.  Additional 


constraints  are  imposed  on  the  density  field  due  to  mass  conservation  and  the  topogra¬ 
phy.  The  resulting  linear  equations  fitting  all  constraints  to  the  data  are  solved  through 
a  minimization  procedure. 

The  horizontal  grid  is  shown  in  Figure  1.  It  is  defined  at  equal  intervals  on  the 
polar  stereographic  projection  used  for  the  domain,  with  a  total  resolution  of  47x25 
grid  points,  and  a  mean  grid  distance  of  ~  127  km.  The  vertical  grid  contains  17 
levels  and  is  stretched,  with  a  resolution  that  varies  from  5m  near  the  surface  to  several 
hundred  meters  near  the  bottom.  The  resolution  is  designed  to  track  variations  in  the 
mixed  layer,  mostly  in  middle  latitudes. 

The  net  flow  (i.e.  barotropic)  model  is  formulated  by  integrating  the  equations  from 
bottom  to  top,  and  taking  its  curl  to  eliminate  the  free  surface  gravity  waves.  The  final 
equation  is  expressed  in  terms  of  the  mass  transport  stream  function  ’P  and  is  as  follows 

^  t?/^(V  •  (1  /H)  •  V)(#  -  *//)  =  V  •  (f0/H)  -  (0/fh2)  ■  (4(*/ff)/dx)  (7) 

+  {9H/ *«).(**/*»)  -  {#H/tiy).{M/#x)  +  V  •  (Ah/H  •  VJ(V*))  (8) 

where  fo  =  f  —  Q,  and  4?  =  ( g/H 2)  •  /  ztrtdz.  Here  f  is  the  surface  stress,  Q  is  the 
bottom  torque,  <xt  is  the  in-situ  density  anomaly,  /3  is  the  NS  variation  of  the  Coriolis 
parameter,  and  Ah  is  the  lateral  viscosity. 

3  Research  Studies  and  Discussions 

We  began  our  modelling  studies  with  the  first  approach  described  above,  in  which  the 
ice  and  the  upper  ocean  are  represented  prognostically  and  the  interior  ocean  diag¬ 
nostically.  The  model  was  initialized  with  temperatures  and  salinities  taken  from  the 
Levitus  climatology,  and  was  forced  by  12-hourly  atmospheric  fluxes  derived  from  the 
NOGAPS  global  atmospheric  prediction  model  of  FNOC  (Fleet  Numerical  Ocean  Cen¬ 
tral  in  Monterey,  CA).  The  integrations  were  performed  for  five  years.  At  the  end  of 
this  period  the  seasonal  behavior  has  converged  to  a  statistical  equilibrium,  such  that 
the  annual  evolutions  of  the  ice  thickness  and  ice  concentration  have  repeated  them¬ 
selves  during  the  fourth  and  fifth  year.  Temperatures,  salinities,  and  velocities  behaved 
similarly.  In  the  fourth  year  the  average  ice  thickness  for  the  month  of  December  has 
evolved  to  the  state  shown  in  Figure  2.  From  the  North  slope  towards  the  North  Pole 
the  ice  thickness  decreases,  as  expected  from  previous  work.  The  NPOC  (Navy  Polar 
Ocean  Center  in  Washington)  analysis  for  ice  edge  location  in  the  Barents  and  Green¬ 
land  Seas  is  indicated  in  the  figure.  The  disagreement  in  the  Barents  Sea  we  attribute 
to  the  weakness  of  the  warm  Norwegian  current,  as  calculated  by  the  diagnostic  interior 
ocean  model  employing  the  (coarse  resolution)  Levitus  climatology.  Not  enough  heat  is 
transported  into  the  ice  edge  region.  The  disagreement  near  Greenland  we  attribute  to 
similar  problems  with  the  Greenland  Current. 

In  figures  3  and  4  we  display  the  temperature  and  salinity  along  the  transect  shown 
in  Figure  1,  during  the  summer  on  day  216.  Three  distinct  regions  can  be  identified: 
under  the  ice,  the  MIZ,  and  the  open  water.  Under  the  ice  the  stratification  tends  to 
be  stronger  than  elsewhere. 

During  the  summer  the  effects  of  the  melting  ice  are  manifested  as  lower  salinities 
under  the  ice  sheet  due  to  salt  dilution.  In  the  winter  months  the  opposite  effect  occurs. 
In  the  MIZ  there  is  warm  water  where  during  the  winter  ice  and  well  mixed  colder  water 
existed.  In  the  open  water  we  see  the  effects  of  the  summer  stratification  of  the  surface 
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layers.  In  the  winter  they  are  well  mixed.  The  salinity  is  high  in  the  open  water  in  the 
summer  and  winter,  reflecting  the  presence  of  warm  saline  Atlantic  water.  In  the  MIZ 
there  is  warm  water  where  during  the  winter  ice  and  well  mixed  colder  water  prevailed. 

In  the  upper  ocean  region,  the  velocity  is  a  superposition  of  the  Ekman,  inertial,  and 
geostrophic  velocities.  An  example  of  it  is  shown  in  figure  5  at  2.5  meters  during  day  8 
of  the  simulation.  The  observed  patterns  are  controlled  by  the  wind  stress  (  transmitted 
through  the  ice  and  applied  to  the  open  water).  The  water  movement  follows  the  ice 
motion  at  an  angle.  The  geostrophic  velocity  varies  from  a  few  cm/sec  to  about  10 
cm/sec.  Its  patterns  reflect  the  Beaufort  gyre,  the  transpolar  drift,  and  the  Greenland 
and  Norwegian  currents.  As  mentioned  before  it  is  calculated  diagnostically  from  the 
Levitus  climotalogy  in  this  approach. 

A's  a  next  step  in  our  development  of  models  of  increasing  complexity,  we  have  intro¬ 
duced  a  prognostic  calculation  of  the  barotropic  ocean  velocity.  Topography  is  included 
in  the  formalism  and  is  coupled  to  the  density  field  via  the  JBAR  representation.  Our 
first  aim  was  to  study  the  evolution  of  the  basic  net  flow  gyres  that  would  arise  without 
coupling  to  the  density  field  or  the  ice.  For  forcing  we  have  used  the  30-year  climato¬ 
logical  winds  of  Dr.  John  Walsh  of  the  University  of  Illinois.  We  have  converged  the 
model  to  a  steady  state  and  studied  the  seasonal  behavior  of  the  gyres  that  evolved.  In 
the  winter  an  anti-cyclonic  gyre  developed  in  the  Beafort  Sea,  together  with  cyclonic 
gyres  in  the  Greenland-Norwegian  and  Barents  Seas.  These  results  are  in  agreement 
with  Semtner(1987)  who  has  also  studied  the  circulation  of  the  Arctic.  In  the  summer 
two  gyres  have  developed  in  the  Arctic  basin,  an  anticyclonic  gyre  in  the  Bering  Sea  side 
and  a  cyclonic  gyre  in  the  Barents  Sea  side  (Figure  6).  This  supports  the  previously 
observed  reversal  of  the  ice  drift  in  the  summer  in  the  Greenland  and  Canadian  sections 
of  the  Arctic  Basin.  It  is  suggested  that  the  climatological  winds  have  a  signature  of 
the  curl  in  them  which  corresponds  to  that  of  the  circulation. 

The  transports  of  the  permanent  gyres. exhibit  a  seasonal  variation.  The  Greenland- 
Norwegian  gyre  has  a  maximum  in  the  winter  and  a  minimum  in  the  summer.  The 
Barents  sea  gyre  has  a  similar  behavior.  The  Beaufort  sea  gyre  has  a  summer  minimum 
with  two  maxima  (Figure  7).  This  seasonal  behavior  agrees  with  Semtner  (1987).  Our 
magnitudes  are,  however,  somewhat  larger  than  his. 
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